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Abstract. For any prime power q and any dimension s 1, a new construction of (t, s)- 
sequences in base q using global function fields is presented. The construction yields an 
analog of Halton sequences for global function fields. It is the first general construction 
of (t, s)-sequences that is not based on the digital method. The construction can also be 
put into the framework of the theory of (u, e, s)-sequences that was recently introduced by 
Tezuka and leads in this way to better discrepancy bounds for the constructed sequences. 



o 

CN ; 1. Introduction 

The construction of low-discrepancy sequences is an important problem in number theory 
and combinatorics and it is also highly relevant for quasi-Monte Carlo methods in scientific 
computing (see [2j|9]). Let us recall the definition of the star discrepancy and the associated 
definition of a low-discrepancy sequence. Let s ^ 1 be a given dimension and let V be a 
point set consisting of the N points x , . . . , Xjv-i in the s-dimensional unit cube [0, 1} S . For 
a subinterval J of [0, 1} S , let A(J; V) be the number of integers n with ^ n ^ iV — 1 and 
x„ G J. Then the star discrepancy D* N (V) of V is defined by 

A(J;T) 



D* N (V)=sup 



N 



X S (J) 



where the supremum is extended over all subintervals J of [0, l] s with one vertex at the origin 
and X s denotes the s-dimensional Lebesgue measure. For a sequence S of points x ,x 1; . . . 
in [0, l] s , we define the star discrepancy D* N (S) for all iV ^ 1 by putting D* N (S) = D* n (Vn), 
! where Vn is the point set consisting of the first iV terms x , . . . , xjv-i of «S. We say that S 
is a low- discrepancy sequence if 

O; (!) D N (S) = 0(N- 1 (\ogN) s ) for all N > 2, 

^ ; where the implied constant is independent of N. This is the smallest possible order of 
magnitude that can currently be obtained for the star discrepancy of a sequence of points in 
[0, l] s . Actually, it is a widely held belief that no smaller order of magnitude can be achieved. 

Historically, the first construction of low-discrepancy sequences for any dimension was that 
of Halton sequences in [3]. For any integer b ^ 2, let Z\, = {0, 1, . . . , b — 1} C Z denote the 



least residue system modulo b. Every integer n ^ has a unique digit expansion 

oo 

n = a r (n)b r 

r=0 

in base b, where a r (n) G for all r ^ and a r (n) = for all sufficiently large r. The 
radical-inverse function (pf, in base b is defined by 

oo 

(j) b {n) = a r (n)b-' r - 1 G [0, 1) for all n G N . 

r=0 

As usual, we write No for the set of nonnegative integers and N for the set of positive integers. 
Now choose pairwise coprime integers b±, . . . ,b s ^ 2 and define the Halton sequence in the 
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bases b\, . . . , b s as the sequence h , hi, . . . with 

K = (<f) bl (n),...,<f) bs (n)) G [0,l) s for n = 0,1,.... 

An explicit upper bound of the form ([T]) for the star discrepancy of this Halton sequence can 
be found in [21 Theorem 3.36]. In practical implementations, the bases bi, . . . , b s are usually 
taken to be distinct prime numbers. 

Modern methods for the construction of s-dimensional low-discrepancy sequences are based 
on the theory of (t, m, s)-nets and (t, s)-sequences which was developed in [7j. This theory has 
a strong combinatorial flavor; see [2j Chapter 6], [HI Chapter 4], and the article on (t,m,s)- 
nets in the handbook pQ. The currently most powerful constructions of (t, s)-sequences use 
global function fields and quite a number of such constructions are available; see [SJ El QUI 
dU EH EH] as well as the surveys in [21 [LS El CE5] • 

In this paper we present a new construction of low-discrepancy sequences using global 
function fields. The construction is comparatively simple and inspired by the construction of 
Halton sequences. It fits very well into the framework of the theory of (u, e, s)-sequences de- 
veloped recently by Tezuka [18J. This theory is an extension of the theory of (t, s)-sequences 
and leads in many cases to better upper bounds on the star discrepancy. The new construc- 
tion is exceptional in the sense that it is the first general construction of (t, s)-sequences 
and (u, e, s)-sequences that is not based on the digital method, which so far was the stan- 
dard method of constructing (t, s)-sequences and (u, e, s)-sequences (see [21 Chapter 4] for a 
description of the digital method). Rather, our construction provides a direct and explicit 
formula for the terms of the sequence (see equation ([7]) in Section [3]). 

The rest of the paper is organized as follows. In Section [2] we present the necessary 
background on (w, e, s)-sequences. Section [3] contains the construction of our Halton-type 
sequences. In Section H] we show that these sequences are both {u, e, s)-sequences and (t, s)- 
sequences with suitable parameters u, e, and t. Upper bounds on the star discrepancy of 
these sequences are discussed in Section [5j 

2. Background on (u, e, s)-sequences 

Tezuka [TS] recently introduced the concepts of (u, m, e, s)-nets and (u, e, s)-sequences. 
For several reasons, e.g., to prove Proposition [1] below and to avoid unnecessary additional 
conditions, his definitions were slightly revised by Hofer and Niederreiter [5]. We follow this 
revised approach. 

Definition 1. Let b ^ 2, s ^ 1, and ^ u ^ m be integers and let e = (ei, . . . , e s ) G N s 
be an s-tuple of positive integers. A (w, m, e, s)-net in base b is a point set V of b m points in 
[0, l) s such that A( J; V) = b m X s (J) for every interval J of the form 

s 

(2) j = n [ ai &-* ( ai +i)&-*) 

i=l 

with integers di ^ 0, ^ a» < b di , and ei\di for 1 ^ i ^ s and with X S (J) ^ 5«- m . 

The classical concept of a (u, m, s)-net in base b corresponds to the special case e = 
(1, . . . , 1) in Definition [TJ For a real number x G [0, 1], let 

oo 

x = Ujb~^ with all yj G Z^, 

3=1 

be a 6-adic expansion of x, where the case yj — b — 1 for all sufficiently large j is allowed. 
For any integer m ^ 1, we define the truncation 

m 

[x}b,m = ^Vjb' 3 - 
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If x = (x^ 1 * , . . . ,x^ s ') G [0, l] s , then the truncation [x]& )m is defined coordinatewise, that is, 

[x]&,m i\p^ ^~\b,mi ■ ■ ■ i \?C \b,m) • 

Definition 2. Let b ^ 2, s ^ 1, and it ^ be integers and let e G N s . A sequence xo, xi, . . . 
of points in [0, l] s is called a (u, e, s) -sequence in base b if for all integers k ^ and m > u 
the points [x n ]b m with kb m ^ n < (k + l)b m form a (it, m, e, s)-net in base b. 

The classical concept of a (it, s)-sequence in base b corresponds to the special case e = 
(1, . . . , 1) in Definition |5J The following result was shown in [3]. 

Proposition 1. Let b ^ 2, s ^ 1, and w ^ 6e integers and let e = (e 1; . . . , e s ) G N s . T/ien 
any (u, e, s) -sequence in base b is a (t, s) -sequence in base b with 

s 

1=1 

For e 7^ (1,...,1), there are currently four families of sequences for which it has been 
verified that they are (it, e, s)-sequences: the Niederreiter sequences constructed in [8], the 
generalized Niederreiter sequences constructed in [17J, the Niederreiter-Xing sequences con- 
structed in pj5], and the Hofer- Niederreiter sequences constructed in [5]. The proof of the 
property of being (it, e, s)-sequences is given in [18] for Niederreiter sequences and general- 
ized Niederreiter sequences and in [5] for Niederreiter-Xing sequences and Hofer-Niederreiter 
sequences. 

3. The construction 

As we have already mentioned, our new construction of low-discrepancy sequences is based 
on global function fields. Recall that a global function field is an algebraic function field 
of one variable over a finite field. We follow the monograph [16] regarding terminology and 
notation for global function fields. 

We write W q for the finite field with q elements, where q is an arbitrary prime power. Let 
F be a global function field with full constant field W q , that is, with ¥ q algebraically closed 
in F. We assume that F has at least one rational place, that is, a place of degree 1. Given a 
dimension s ^ 1, we choose s + l distinct places P^, Pi, . . . , P s of F with deg(Poo) = 1. The 
degrees of the places Pi, . . . , P s are arbitrary and we put e; = deg(p) for 1 ^ i ^ s. Denote 
by Op the holomorphy ring given by 

o F = p| Op, 

where the intersection is extended over all places P ^ P^ of F and Op is the valuation ring 
of P. For each % = 1, . . . , s, let jpj be the maximal ideal of Op corresponding to Pj. Then the 
residue class field Fp t := Of/ pi has order q £t (see [16, Proposition 3.2.9]). We fix a bijection 

(3) a Pi : F Pi -> Z qH . 

For any divisor D of F , we write C(D) for the Riemann-Roch space associated to D and 
£(D) for the dimension of the vector space C(D) over W q (see [161 Section 1.4]). For each 
i — 1, . . . , s, we can obtain a local parameter Zi G Op at pi, by applying the Riemann-Roch 
theorem [161 Theorem 1.5.17] and choosing 

Zi G ^(fcPoo - P) \ CikP^ - 2P) 

for a suitably large integer k. Then for every / G Of, we have a local expansion of / at pi 
of the form 

oo 

(4) f = Y,M with all f id eF Pi . 

3=0 
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Compare with [T5J pp. 5-6] for a detailed description of how to obtain a local expansion. 
Now we define the map : Op —> [0, l] s by 

(oo oo \ 

j=0 j=0 J 

The map <fi can be viewed as an s-dimensional radical- inverse function for Of- 
We arrange the elements of Op into a sequence by using the fact that 

oo 

Of = |J ^(mPoo). 

m=0 

The terms of this sequence are denoted by / , ft, . . . and they are obtained as follows. Con- 
sider the chain 

C(0) C £(P 00 ) C £(2P 00 ) C • • • 

of vector spaces over ¥ q . At each step of this chain, the dimension either remains the same 
or increases by 1. From a certain point on, the dimension always increases by 1 according to 
the Riemann-Roch theorem. Thus, we can construct a sequence Vq,v±, . . . of elements of Of 
such that for each m G No, 

{«0,Vlj • • • ,Vi( mPoa )-l} 

is an F 9 -basis of £(raPoo)- For n G N , let 



yja r (n)(f with all a r (n) G Z 9 



r=0 

be the digit expansion of n in base q. Note that a r (ri) = for all sufficiently large r. We fix 
a bijection 77 : Z 9 — > ¥ q with 77(0) = 0. Then we define 

(6) f n = r)(a r (n))v r G Op for n = 0, 1, . . . . 

r=0 

Note that the sum above is finite since for each n G No we have r)(a r (n)) = for all sufficiently 
large r. 

Using (j^J) and (jBJ), we now define the sequence x , Xi, . . . of points in [0, l] s by 

(7) x n = 4>(f n ) for n = 0, 1, . . . . 

This is our Halton-type sequence obtained from the global function field F. We will show 
that this sequence is a (it, e, s)-sequence in base q for some u and e (see Theorem [TJ below) 
and also a (t, s)-sequence in base q for some t (see Corollary [T] below). 

Remark 1. The sequence x , x 1; . . . is obtained by the direct and explicit formula (j7j) and it is 
easy to see that, in general, it cannot be produced by the digital method for the construction 
of (t, s)- sequences. The simple reason is that the digital method for the construction of (t, s)- 
sequences in base q is a linear-algebra technique which works entirely in the finite field ¥ q 
(compare again with [21 Chapter 4]), whereas our construction uses extension fields of ¥ q 
of degrees ei,...,e s . Moreover, whenever ei ^ 2 there is no reasonable way in which the 
bijection ap. in ([3]) can be interpreted as an F^-linear map. There is only one situation in 
which our construction can be put into the framework of the digital method, namely when 
q is a prime number and = 1 for 1 ^ i ^ s. In this case the bijections a Pi in Q can be 
chosen as identity maps. 

Remark 2. Tezuka [T7j introduced a construction of Halton-type sequences based on rational 
function fields over finite fields. If we specialize our construction to the case where F is the 
rational function field over ¥ q , then our construction is in general different from Tezuka's 
construction. One significant difference is that Tezuka's construction can be described in 
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terms of the digital method, whereas our construction cannot in general be put into the 
framework of the digital method (see Remark [T]). 

4. The main result 

Before we formulate our main result in Theorem [1] below, we prove the following lemma. 
We keep the notation from the previous sections. 

Lemma 1. Let g be the genus of F and let k ^ and m ^ g be integers. Let the elements 
f n G Op be defined by (jBj). Then for f G F we have f = f n for some integer n with 
kq m ^ n < (k + l)q m if and only if f = h + c with h G Of depending only on k and m and 
with c G C((m + g — l)Poo). 

Proof. The integers n with kq m ^ n < (k + l)q m have digit expansions in base q of the form 

m—1 oo 

n = a r (n)q r + a r (n)q r , 

r=0 r=m 

where the a r (n) with ^ r ^ m — 1 range independently over Z q and the a r (n) with r ^ m 
depend only on Thus, the corresponding f n have the form 

m—1 oo m—1 

(8) f n =y j n{a r {n))v r + ^2 v{ a r{n))v r =: 2J a r v r + 

r=0 r=m r=0 

where the a r = rj(a r (n)) with ^ r ^ m — 1 range independently over ¥ q and h £ Op 
depends only on k and m. 

Now we consider C((m + g — 1)-Poo) and note that deg((m + g — l)Poo) = m + g — 
1 ^ 2g — 1. Therefore £((m + g — l)Poo) = m by the Riemann-Roch theorem. It follows 
that {vo,vi, . . . ,f m _i} is an F 9 -basis of £((m + g — 1)P 00 ), and so the linear combinations 
J^^ 1 ot r v r range exactly over £((m + g — l)Poo). In view of (jSJ), the proof is complete. □ 

Theorem 1. The sequence x , xi, . . . given by is a (u, e, s) -sequence in base q with u = g 
and e = (ei, . . . , e s ), where g is the genus of F and = deg(p) for 1 ^ i ^ s. 

Proof. We proceed by Definition [2] and fix integers k ^ and m > g. We have to show that 
the point set Vk, m consisting of the points [x n ] ?jm = [<f>(f n )]q,m with kq m ^ n < (k + l)q m 
forms a (g, m, e, s)-net in base q. 

Let J be an interval of the form ([2]) with A S (J) ^ q 9 ~ m , that is, 

s 

(9) £di<m-0. 

i=i 

Note that we assume ej|dj for 1 ^ i ^ s according to Definition [TJ Let / G with local 
expansions as in (J3J). Then [<j)(f)} qtm G J if and only if 

oo 

E^(/m-i)(^)- J 
This is equivalent to 

di/ei 

Wi := °r(.i)j i)'/ " e Kg" di , (flj + i)g~ di ) forU^s. 
i=i 

This condition is satisfied if and only if, for 1 ^ i ^ s, the first di g-adic digits of Wi agree 
with the first di g-adic digits of aiq~ di and thus have fixed values depending only on J. Note 



G [aiq~ d \ [at + l)q~ di ) for 1 < i < s. 
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that these d{ fixed g-adic digits of u>, correspond to rfj/e, fixed g £i -adic digits. Therefore for 
some dij G Z q H, 1 ^ J ^ d%/ e h 1 ^ 2 ^ s, we get the condition 

<7Pi{fi,j-i) = <Hj for l^j^di/d, 1 ^ i ^ s, 
or the equivalent condition 

(10) jij-x = &p?{a>i,j) =■ G F Pt for 1 ^ j ^ di/e i; 1 ^ i ^ s. 
Now we recall 01]) and put 

(di/ei)-l 

/ij = f3ij + iz{ G Op for 1 ^ i ^ s. 

j=0 

Then f jT0|) is equivalent to 

(11) / = ^ mod pf /ei for 1 ^ z < s. 

Note that the elements h\, . . . , h s depend only on J. 

We consider now the number A(J; 7\ m ) of integers n with kq m ^ n < (k + l)q m such that 
[<fi(fn)]q,m ^ J ■ As we have just seen, the latter condition is equivalent to ffTTj) holding for 
/ = f n - By Lemma [1] the condition kq m ^ n < (k + l)q m is equivalent to f n = h + c with 
fixed /i G Op and some c G £((m + g — l)Px>)- Thus, both conditions hold simultaneously if 
and only if c G £((m + g — l)Poo) satisfies 

(12) c=hi-h mod p* /e * for 1 ^ i ^ s. 
Consider the map 

^ : C ((m + 5 - 1)^) ^R:= (0 F /p? l/ei ) x • • • x {Op 

defined by 

ip(c) = (c mod pf l/e \ . . . , c mod pf /es ) for c G £((m + # - l)Poo). 

Note that ip is a linear transformation between vector spaces over ¥ q . We claim that ip is 
surjective. To prove this, it suffices to show that 

dim(£((m + g - ljP^/ker^)) = dim(P). 

The dimension dim(P) is easily computed to be 

s s 

dim(P) = ^dim(0 F /pf /ei ) =J2 d i- 

i=l i=l 

From the definition of r/> it is clear that 

s 

(13) ker(V) = £((m + 5 - 1)? M - J^/e^p). 

i=i 

Furthermore, by taking into account we get 

s s 

(14) deg ((m + <? - VjP^ - J^/e^p) = m + <?-l-J] C f i ^2<7-l. 

j=l i=l 

Therefore 

s 

dim(£((m + g- l^/ker^)) = dim (c((m + g- l)P x )/c(fm + g- l)Poo - J^/e^P 

= deg(^(d i /e,)P i ) =^di 

8=1 4=1 

by the Riemann-Roch theorem, and so ■?/> is indeed surjective. 



i=i 
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Since ip is surjective, the number of c G £((m + g — 1)Pqo) satisfying the system of con- 
gruences in (TT2"j) is equal to #ker (?/>). By (TIB"]) . fTTTl) . and the Riemann-Roch theorem we 
get 

s 

dim(ker(V)) = l({m + g - l)Poo - J^/ e i) P *) 

i=i 

= deg f(m + #- l)Poo - ^(di/e^PiJ + 1- # = to-^cZ;. 

j=i i=i 

Thus, we have shown that 

A(J;V k , m ) = #ker(^) = q™~ H=i* 

On the other hand, we have 

g m A s (J) = g m -^=i* 

by the form of J, and so 7\ m is indeed a (g, m, e, s)-net in base q according to Definition[TJ □ 
Corollary 1. The sequence x ,X!, . . . given by is a (t, s)- sequence in base q with 

s 

t = g + Y,( e i- 1 h 

i=l 

where g is the genus of F and e, = deg(Pj) for 1 ^ i ^ s. 

Proof. This follows from Theorem [JJ and Proposition [TJ □ 

Remark 3. The Niederreiter-Xing sequences constructed in [19] use the same ingredients as 
our construction, to wit a global function field F with full constant field ¥ q and s + 1 distinct 
places Poo, Pi, . . . , P s of F with deg(Poo) = 1 and the degrees e, = deg(p) for 1 ^ i ^ s 
being arbitrary. It was shown in [J_H] that a Niederreiter-Xing sequence with these data is a 
(t, s)-sequence in base q with the same value of t as in Corollary [TJ Furthermore, Hofer and 
Niederreiter [5] proved that the same Niederreiter-Xing sequence is also a (g, e, s)-sequence 
in base q with e = (ei, . . . , e s ), that is, with the same parameters as in Theorem [TJ Thus, the 
Niederreiter-Xing sequence and our new Halton-type sequence have a similar distribution 
behavior. However, it should be pointed out that our new construction is considerably 
simpler than the one in [19]. Since the Niederreiter-Xing sequences and our new Halton- 
type sequences have the same parameters, the usual procedure of optimizing the t-value for 
Niederreiter-Xing sequences (see [12] and [13, Section 8.3]) applies in the same way to our 
new Halton-type sequences. Consequently, our new construction yields (t, s)-sequences in 
base q such that, for fixed q, the parameter t grows linearly in s as s — > oo and is therefore 
asymptotically optimal. 



5. Discrepancy bounds 

The star discrepancy of any (t, s)-sequence S in any base 6^2 satisfies an upper 

bound of the form (TTJ). In fact, we have 

(15) D* N {S) ^ CN~\logN) s + OiN^ilogN) 3 - 1 ) for all N ^ 2 

with the constant C > and the implied constant in the Landau symbol depending only on 
b, s, and t. This was first shown in [TJ Section 4] (see also [HI Theorem 4.17] for a convenient 
formulation). The currently best values of C are those of Faure and Kritzer [3J, namely 
C = C*fk given by 

c = (-S) if & is even, 

if 6 is odd. 



cS 
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(16) C FK ^ 



It follows that for every base b ^ 2 we have 

b l 1/6-1 
s! ' 2 V2log6 / 

Tezuka [18] established an upper bound on the star discrepancy D* N (S) of any (u,e,s)- 
sequence S in base b. This bound is also of the form ( ITS]) with the constant C = Cxez given 
by 

6" * L6 e V2j 6" fe ei+ - +es 
1 J Tcz si Al ei l 0g 6 ^ s! (21og6) s ei ---e s ' 

Note that our Definition [2] of a (w, e, s)-sequence in base b is slightly stronger than Tezuka's 
original definition in [18] , and so Tezuka's discrepancy bound is obviously valid for our concept 
of a (u, e, s)-sequence in base b as well. 

Our new Halton-type sequences are (t, s)-sequences in base q as well as (u, e, s)-sequences 
in base q with suitable t, u, and e. Therefore we can apply two versions of the discrepancy 
bound ( 1151) . namely the one with C = Cyk an d the one with C = C^ cz . It is of interest to 
compare these two bounds. Recall from Theorem [T] and Corollary [T] that for these sequences 
we have u = g, the genus of the global function field F, as well as e = (ei, . . . ,e 8 ) and 
t = g + J2i=i( e i ~ I); where = deg(Pj) for 1 ^ i ^ s. Then from (TT6]) and (fT7j) we obtain 

i=i 



C Tez 2 V g 

Thus, if 



n 



e { > 2 
1 



then the discrepancy bound with the constant Ct cz is better than the one with the constant 
Cyk- The condition f TT8]) will be satisfied in many cases. For instance, it was shown in [5] 
that if we arrange all distinct places 7^ of F into a list Pi, P2, . . . in an arbitrary manner, 
then for any 0<£<l/q , we have 



n 



»((l/9)-«0« 



deg(P) > (log,s) 

i=l 

for all sufficiently large s, where log 5 denotes the logarithm to the base q. Hence in this 
situation the condition f lTSl) is satisfied for all sufficiently large s, and so the discrepancy 
bound ffTBT) with the constant C = C^ cz is better than the one with the constant C = Cfk- 

The discussion above shows that the theory of (u, e, s)-sequences is a useful extension of 
the theory of (t, s)-sequences since, for instance, it can lead to improved discrepancy bounds 
in comparison with the discrepancy bounds for (t, s)-sequences. 
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